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We identify a route towards achieving a negative index of refraction at optical frequencies based 
on coupling between plasmonic waveguides that support backwards waves. We show how modal 
symmetry can be exploited in metal-dielectric waveguide pairs to achieve negative refraction of both 
phase and energy. By properly controlling coupling between adjacent waveguides, a metamaterial 
consisting of a one-dimensional multilayer stack exhibiting an isotropic index of —1 can be achieved 
at a free-space wavelength of 400 nm. The general concepts developed here may inspire new low-loss 
metamaterial designs operating close to the metal plasma frequency. 
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Metamaterials allow control over the propagation of 
electromagnetic waves in ways beyond those provided by 
naturally occurring materials. Since the prediction that a 
material with a negative index of refraction can be used 
to construct a perfect lens pQ, and spurred on by the 
paradigm of transformation optics [2 J, a quest for me- 
dia that possess a negative index has resulted in a wide 
range of metamaterials operating in different frequency 
regimes [3]. Many designs of such 'left-handed' materials 
are based on subwavelength metallic resonant scattering 
elements that exhibit a simultaneously negative electric 
and magnetic response to light. These geometries can 
be related to archetypical structures such as split rings 
at low frequencies, where e m <C — €d (s m and e^ are the 
dielectric constants of the metal and its dielectric sur- 
rounding, respectively) |3]. However, using this route 
it has proven to be very difficult to achieve a negative 
index with acceptable losses at visible frequencies, i.e., 
close to the metal plasma frequency uj p [U{6]. Moreover, 
most designs only show the desired behavior for a lim- 
ited range of angles. These limitations are less severe in 
indefinite materials, which rely on hyperbolic dispersion 
to achieve negative refraction of energy [TUTU] . Because 
the phase index in these materials is positive, indefinite 
media however do not possess all the same properties as 
left-handed media. In a left-handed medium the direc- 
tion of energy flow is opposite to that of the wavevector 
(k • S < 0) [11] [12] . An ideal negative index metamate- 
rial has a phase index that is independent of angle, which 
results in the wavevector k and Poynting vector S being 
entirely antiparallel. Both phase and energy are then 
refracted negatively at an interface of a positive index 
medium with the metamaterial. 

In this work, we propose a strategy to design meta- 
materials that have a three-dimensional negative index 
of refraction at visible frequencies. The working princi- 
ple relies on coupling between plasmonic waveguides that 
are arranged in a stacked geometry. We use the optical 
properties of metals close to the plasma frequency, where 



— £d < £m < 0, which allow the constituting waveguides 
to support backwards surface plasmon polariton (SPP) 
waves [T3HT8] . We show that the symmetry of the waveg- 
uide modes is crucial to the realization of negative refrac- 
tion in the metamaterial. By exploiting symmetric modes 
in arrays of plasmonic waveguide pairs, the coupling be- 
tween adjacent waveguides can be tailored to achieve 
negative refraction of both wavevector and energy. The 
design strategy presented in this manuscript allows the 
attainment of an isotropic negative index of refraction 
of —1 for TM-polarized light in a one-dimensional array 
of plasmonic slab waveguides. The material operates at 
visible frequencies with comparatively small absorption 
losses. 

Wave propagation in an array of coupled waveguides 
can be described by coupled mode theory in the weakly 
coupled paraxial limit [9j [T9] , We consider a material 
comprised of slab waveguides aligned along the x?/-plane 
and stacked in the z direction (as sketched in the inset 
of Fig. [I]). In the xy-plane, parallel to the slabs, the 
wavevector is independent of angle and its magnitude j3 r 
will be slightly altered with respect to the mode wavevec- 
tor f3 in a single isolated waveguide due to the presence of 
the adjacent waveguides. The angle at which a beam is 
refracted that is incident in the xy-plane from a homoge- 
neous medium onto an interface parallel to the ?/z-plane 
can therefore be found from Snell's law, using the refrac- 
tive index /3'c/lj. For the material to be left-handed and 
refraction of both wavevector and energy to be negative 
in the x?/-plane, the mode wavevector /3 of the constituent 
waveguides must therefore be negative. 

In the xz-plane, propagation and refraction is governed 
by waveguide coupling. Taking only nearest-neighbor 
interactions into account, the amplitudes in adjacent 
waveguides are related through [T9] 



dcij/dx = ipctj + in (a 
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where clj is the complex amplitude of the mode in the 
j th waveguide and k a coupling constant. Solving Eq. (Ly) 
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FIG. 1: (color) Wavevector diagram showing isofrequency 
contours of a homogeneous positive index medium (blue) and 
left-handed waveguide array with k > (red) or k < (pur- 
pie). 



yields a simple relationship between the wavevector com- 
ponents parallel (k x ) and normal (k z ) to the waveguides: 



k x = f3 + 2k cos (k z d) , 



(2) 



where d is the center-to-center distance between adjacent 
waveguides. This relationship is sketched in a wavevec- 
tor diagram in Fig. [I] for different signs of the coupling 
constant k and taking /3 < 0. Also indicated is the 
isofrequency contour for a homogeneous positive index 
medium (blue curve). Since the problem of interest is 
that of refraction of a plane wave incident in a homoge- 
neous medium from x = — oo (see inset of Fig. IT]), only 
curves are drawn for which S x > 0. By definition, the 
group velocity (and in the absence of loss also S) is ori- 
ented normal to an isofrequency contour. Given a par- 
ticular incident wavevector k 1 , the refracted wavevector 
k r can be derived from the wavevector diagram by con- 
serving the wavevector component k z (up to an integer 
times 2ir/d). We can see from Fig. [I] that k r • S r < 0, 
which results from the fact that (3 < 0. This material can 
therefore rightly be called left-handed, and wavefronts 
will thus be refracted negatively in the xz-plane as well 
as in the xy-plane. Importantly, however, whether the 
refraction of energy, given by the direction of S r , is posi- 
tive or negative is determined by the sign of the coupling 
constant n. If k, < 0, refraction of energy is negative 
(i.e., SI SI < 0), whereas it is positive if n > 0. Three- 
dimensional negative refraction of both phase and energy 
can therefore only be realized when both (3 and k, are neg- 
ative. 

Guided modes with /3 < can be called 'backwards' 
since energy and phase are counterpropagating [TT] . 
They can exist in plasmonic waveguides in which more 
energy is guided in metal than in dielectric, since inside 
the metal S is oriented opposite to that in the dielec- 
tric. The antisymmetric SPP mode guided by a thin 
dielectric slab clad with metal, termed a metal-dielectric- 
metal (MDM) waveguide, can acquire a backwards char- 
acter under certain conditions, as pointed out by sev- 
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FIG. 2: (color) Dispersion in MDM waveguides for varying 
dielectric thickness t, showing the appearance of backwards 
SPPs for small t. On the right, characteristic E z profiles 
of the symmetric (s) and antisymmetric (a) SPP modes are 
displayed. 



eral authors [T3HT8] . Figure [2] depicts the dispersion 
curves of the propagating symmetric and antisymmet- 
ric SPP modes in a MDM waveguide for various dielec- 
tric thicknesses t, calculated as described in [18 using 
£m — 1 ~ UJ p/ uj2 an d £d = 2.25. For simplicity, we neglect 
absorption here, but as we will see the same physics ap- 
plies once it is included. The symmetry of the transverse 
electric field assigns the mode symmetry. The antisym- 
metric mode acquires a negative group velocity du/d/3 
for thicknesses t < (ttc/2uj p ) y/(l + £<*) /s&, when its cut- 
off frequency uj c = (2c j ' y/e&t) tan -1 yj — e m /ed becomes 
larger than the surface plasmon resonance frequency 



^sp 



LJp/y/1 + £d- I n the absence of absorption, this 
directly means that energy and phase are counterprop- 
agating. Because for uj sp < uj < uj p the antisymmetric 
mode is the only propagating mode, the waveguide can 
be described by the effective index /3c/ lj, where we have 
to use the 'mirrored' branch with j3 < when we assume 
energy propagates forwards [18] . This can lead to nega- 
tive refraction of light into such a waveguide [16, 18)120] . 
This effect is however purely two-dimensional, restricted 
to the plane of the waveguide. 

Now we must ensure that these modes, when coupled 
in a three-dimensional array, lead to a negative coupling 
constant k, to allow negative refraction of energy. The 
coupling constant can be approximated from the modal 
fields in an isolated waveguide by [21 



k = 



_ (e d -e m )J d E^(z)-E(z-d)dz 
U;c ° r °° x • (Et x H + E x Ht) dz ' 
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The use of the adjoint fields E"^ and H"^ (obtained by 
substituting —j3 for f3 and —uo for uo) is important for 
the correct treatment of absorbing media [2TJ [22] . The 
integral in the numerator of Eq. (|3| is performed over 
the dielectric region of the waveguide only [22] . The sign 
of this mode overlap integral determines the sign of n. 
Considering the mode profile E z (z) and the displaced 
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FIG. 3: (color) (a) E z (z) profile of the antisymmetric mode in 
a MDM waveguide (red), together with E z (z — d) (blue), (b) 
Refraction into an array of MDM waveguides. Arrows depict 
the direction of energy flow, (c) E z profile of the backwards 
symmetric mode in a pair of MDM waveguides, (d) Refraction 
into an array of strongly coupled MDM waveguide pairs. 



field E z (z — d) (shown in Fig. [3ja)) of the antisymmet- 
ric mode in a MDM waveguide, we see that the overlap 
integral J\ E\ (z) E z (z — d) dz is positive. This is the re- 
sult of the fact that (1) E z switches sign across a metal- 
dielectric interface and (2) the mode is antisymmetric. 
A similar reasoning proves that J d E]. (z) E x (z — d) dz is 
positive as well [22 . As a result, a stack of equally spaced 
MDM waveguides as depicted in Fig.[3Fb) will in fact have 
n > and therefore exhibit positive refraction of energy 
in the xz-plane. It therefore does not constitute an ideal 
negative index material allowing negative refraction of 
both phase and energy in all directions. 

As outlined above, the reason that k, is positive is re- 
lated to the antisymmetric nature of the mode in a MDM 
waveguide. To achieve negative refraction, a backwards 
mode with a symmetric nature is needed instead. This 
exists for example in a pair of strongly coupled MDM 
waveguides, in which eigenmodes resemble even and odd 
superpositions of the antisymmetric modes in the indi- 
vidual waveguides. The odd superposition has a mode 
wavevector (3 that is slightly smaller than that of a single 
waveguide and also negative [23] . Most importantly, this 
mode has a symmetric field profile (as shown in Fig.[3|c)), 
which causes the coupling constant n in Eq. O) to be neg- 
ative for coupling between adjacent waveguide pairs. An 
array of weakly coupled pairs of strongly coupled MDM 
waveguides, i.e., a metal-dielectric multilayer stack with 
alternating thick and thin metal layers as depicted in 
Fig. [3Fd), will therefore exhibit negative refraction of 
both energy and wavevector in the xz-plane. An im- 
portant second advantage to using the symmetric mode 
in a pair of MDM waveguides is that it can be excited 
by light incident along the x axis from outside the meta- 
material, whereas the antisymmetric mode has vanishing 
mode overlap with such a wave. 
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FIG. 4: (color) (a) Calculated wavevector diagrams for vary- 
ing metal thickness between waveguide pairs, for the param- 
eters described in the text. The red dotted curve shows 
the isofrequency contour for a homogeneous medium with 
n = ±1. (b) Three-dimensional wavevector diagram show- 
ing isofrequency surfaces for three different values of Ao, to 
show that n is isotropic and negative, (c) Snapshot of the H y 
field distribution showing refraction of light from air into the 
metamaterial extending to the right of the vertical line, as cal- 
culated with FDTD. The arrows show the unit-cell-averaged 
direction of the Poynting vector. 



An ideal negative index metamaterial has an isotropic 
phase index, which results in a circular isofrequency con- 
tour. The curvature of the isofrequency contour in Fig. [I] 
is related to the coupling strength k, which can be con- 
trolled by varying the thickness t m of the metal layer sep- 
arating adjacent pairs of waveguides. Figure [4^ a) shows 
isofrequency contours for TM polarization for varying t m , 
calculated using a transfer matrix formalism. The free- 
space wavelength Ao is 400 nm. We use the real part of 
the dielectric constant of Ag (e m = —4.43 [24]). We take 
^fsa = 3.2, comparable to TiC>2 at this wavelength. A 
large dielectric index ensures uj sp is relatively small, ex- 
tending the operating frequency window to a large part of 
the visible spectrum. The thicknesses of both dielectric 
layers and of the thin Ag film are 28 and 35 nm, re- 
spectively. Reducing t m increases the coupling strength, 
resulting in a stronger modulation of k x . Interestingly, 
the isofrequency contours evolve from the cosinusoidal 
behavior described in the weak coupling limit by Eq. [2] 
to almost fully circular for t m = 75 nm. Since the struc- 
ture is isotropic in the £?/-plane, the isofrequency surface 
plotted in a three-dimensional wavevector diagram now 
resembles a sphere, as shown in Fig. plb). Wave prop- 
agation in the metamaterial can thus be described by 
an isotropic negative index of refraction n that deter- 



mines the refraction angle of both phase and group ve- 
locity. The magnitude of n at a given frequency can be 
controlled most effectively by changing j3 through vary- 
ing the dielectric thickness (see Fig. |2j. In this case, 
it has been chosen to result in n = — 1. In Fig. plb), 
we have also indicated the isofrequency surfaces for two 
slightly different frequencies. The wavevector magnitude 
decreases for increasing frequency, confirming that the 
Poynting vector is oriented inwards and the material is 
left-handed. 

It is interesting to note that for strong enough cou- 
pling, even propagation normal to the interfaces (along 
z) is allowed. The geometry and physical mechanism 
is different than for metal-dielectric multilayer stacks in 
the hyperbolic dispersion regime [25]. The existence of 
transparency windows in metal-dielectric multilayers is 
well known, and can be attributed to resonantly coupled 
Fabry-Perot cavities [26]. In the currently proposed sys- 
tem, the fundamental harmonic of the Bloch wave prop- 
agating along z can be described by a negative index. Of 
course, the fact that the wave comprises many Bloch har- 
monics has implications for the application of this meta- 
material. For example, when used to construct a perfect 
lens, the largest magnitude of the wavevector component 
k z that can be transferred by the lens is ir/d, i.e., the 
resolution along z is limited by the dimension of a single 
waveguide. This can however be designed to be signifi- 
cantly smaller than Ao/2. 

When introducing realistic absorption losses in the pre- 
vious example (e m = —4.43 + 0.2H [24]), Re (n) is not 
notably affected [22 J. The figure of merit, defined as 
Re(|k|) /Im(|k|), varies from 14.7 along x to 8.7 along 
z. Considering the fact that the operating wavelength is 
only 400 nm in this metallodielectric metamaterial, ab- 
sorption is thus comparatively small [5, 6, 27 , and not 
strongly dependent on the propagation direction. 

To illustrate negative refraction of light into the meta- 
material, Fig. [4Fc) shows the result of a finite- difference 
time-domain simulation in which the interface of the 
metamaterial truncated along a plane parallel to yz is 
illuminated with a TM-polarized plane wave incident un- 
der a 30° angle from air. The material parameters are 
the same as before, and include absorption in Ag by lo- 
cally fitting the dielectric constant to a Drude model 
with dissipation. The arrows indicate the unit-cell av- 
eraged direction of S, showing that not only the wave- 
fronts, but indeed also the energy is refracted negatively, 
with n = — 1. The fraction of energy coupled to the re- 
fracted wave in the metamaterial is 2%. However, in this 
example impedance matching was not optimized. Simu- 
lations have shown that coupling efficiencies up to 20% 
are readily available by tuning the layer thicknesses [22] . 
Moreover, structuring the interface may further improve 
impedance matching. Further studies should identify the 
limits and importance of impedance matching in a prac- 
tical application, as well as derive values of effective per- 



mittivity and permeability for this metamaterial. 

In conclusion, we have identified a route towards the 
construction of (left-handed) negative index metamateri- 
als using plasmonic waveguide coupling. We underlined 
the role of backwards SPP modes and modal symmetry, 
and showed that the coupling can be controlled to achieve 
negative refraction of both energy and wavevector. A 
carefully designed one-dimensional stack of metallic and 
dielectric layers with can be characterized by a single 
three-dimensional isotropic negative index of refraction 
of n = —1 at a wavelength of 400 nm. In this respect, 
it differs strongly from photonic crystal slabs exhibiting 
negative index behavior [28] [29], which are essentially 
two-dimensional and in which backward wave behavior 
relies fully on Bragg diffraction. The structure proposed 
here only transmits TM-polarized light, but polarization- 
independent behavior could be found in designs that ex- 
ploit alternative waveguide geometries p~4j [3D]. Using 
Ag, the operating frequency of these metamaterials ex- 
tends from the visible to the near-UV regime, with fig- 
ures of merit of the order of 10. The same mechanisms 
can be applied throughout the electromagnetic spectrum, 
using conductive oxides in the near-infrared [31], SiC in 
the mid- infrared [14] , (doped) semiconductors at infrared 
and THz frequencies [51] , and wire- mesh metamaterials 
in the microwave regime [52] . We envisage that the con- 
ceptual framework outlined here may lead to a new class 
of metamaterial designs exploiting the properties of met- 
als near the plasma frequency. 
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SUPPLEMENTARY INFORMATION 

Modal symmetry and the sign of the coupling 
constant 

As argued in the main text, whether energy is refracted 
positively or negatively is controlled by the sign of the 
coupling constant k and the mode wavevector /3. The 
coupling constant is approximated in the weak coupling 
limit by Eq. (3). A more general expression for the cou- 
pling strength between a waveguide aligned along x with 
dielectric constant profile e (z) and an identical adjacent 
waveguide spaced at a distance d is pQ 



: UJSQ- 



f^ oo [e(z)-e(z-d)]EHz)^E(z-d)dz 



/^x-(Et xH + ExHt)^ 



(SI) 

where E and H are the mode fields of a single waveg- 
uide. The use of the adjoint field instead of the com- 
plex conjugate allows the dielectric constants to be com- 
plex [1 j. The fact that the integrand in the numerator 
scales with [e (z) — e (z — d)] can be understood when re- 
lating the coupling constant to the perturbation of the 
field in one waveguide by the presence of the adjacent 
waveguide, with the difference of the dielectric constant 
profiles taking the role of a potential difference. In the 
MDM waveguides considered in this work, this difference 
vanishes everywhere except in the dielectric regions, and 



Eq. (SI) reduces to Eq. (3). Apart from brute- force eval- 
uation of Eq. (3), the sign of n can also be found from 
simple arguments. In the case of lossless media, the ad- 
joint field E^ could be replaced by the complex conju- 
gate E* and the denominator of Eq. (3) then equals four 
times the total power carried by the waveguide, which 
is evidently a positive quantity. The factor (e<± — e m ) is 
positive as well. Therefore, the sign of k, is determined by 
the overlap integral of the modal electric field E (z) with 
the electric field in a neighboring waveguide E(z — d), 
evaluated over the dielectric region. 

As an example, we first consider the symmetric mode 
in a MDM waveguide at frequencies uo < cj sp , for which 
both /? and v g are positive. As sketched in Fig. [5] (Sup- 
plemental) (a), the overlap of the E z (z) component in the 
dielectric core with the field E z (z — d) in the cladding of 
the adjacent waveguide (j d E\(z) E z (z — d) dz) is neg- 
ative, since E z switches sign across a metal-dielectric 
boundary. We used the symmetry of both E z and H y 
to label this mode as being 'symmetric'. In contrast, the 
z profile of the longitudinal component E x of a symmet- 
ric mode is itself antisymmetric (and vice versa). How- 
ever, unlike E z , E x does not change sign across a metal- 
dielectric interface, and hence the net overlap integral 
f d El (z) E x (z — d) dz < 0. This can be judged from the 
mode profile of E x in Fig. [5] (Supplemental) (b). Note 
that we plot Im (E x ) since E x (x, z) is imaginary for the 
particular value of x that was used to calculate the (real) 
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FIG. 5: (Supplemental). Coupling and refraction of the sym- 
metric mode in the positive-index regime, (a) Typical E z (z) 
mode profile of the symmetric mode in a MDM waveguide 
for cj < cj sp (red), and the displaced field E z (z — d) (blue), 
showing how the mode overlap in the dielectric core is nega- 
tive, (b) E x (z) (red) and E x (z — d) (blue) mode profile, with 
also negative overlap, (c) Sketch of the wavevector diagram 
with the isofrequency contour of a homogeneous medium from 
which a plane wave with wavevector k 1 is incident (blue) , and 
the isofrequency contour of the symmetric mode in a coupled 
stack of MDM waveguides, with (3 > and k, > (red), (d) 
Refraction of energy is negative in the xz-plane, but positive 
in the xy-pl&ne. 



E z profile in Fig. [5] (Supplemental) (a). Since the contri- 
butions of E z and E x to the overlap integral in Eq. (3) 
are both negative, n is always negative. As a result, the 
isofrequency contour (Eq. (2)) in the xz-plane of a stack 
of MDM waveguides for uo < uo sp is such that energy is re- 
fracted negatively [2 , as can be seen from the wavevector 
diagram in Fig. [5] (Supplemental) (c). The isofrequency 
contour resembles hyperbolic dispersion for small k z , and 
the arguments above could be similarly applied to explain 
negative refraction in other indefinite media [3]. How- 
ever, refraction of the wavevector is clearly positive and 
angle-dependent in this right-handed metamaterial, and 
refraction of both phase and energy in the £?/-plane is 
exclusively positive as well, as sketched in Fig. [5] (Sup- 
plemental) (d) [4]. 

The antisymmetric backwards mode in a thin MDM 
waveguide for uo > uo sp exhibits opposite behavior com- 
pared to the symmetric mode in many respects. The 
field E z (z — d) in the cladding of the displaced waveg- 
uide now has a net positive overlap with the field E z (z) 
in the core of the waveguide due to the antisymmetric na- 
ture of this field component, as was clear from Fig. 2(b). 
Of course, the longitudinal component E x should be an- 
alyzed in a similar fashion. In fact, for the backwards 
waves considered here the longitudinal component gen- 
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FIG. 6: (Supplemental), (a) Field profile of longitudinal elec- 
tric field E x (z) of the antisymmetric mode in a MDM waveg- 
uide (red), as well as the displaced field E x (z — d) (blue), 
having positive overlap, (b) Profiles of E x (z) (red) and 
E x (z — d) (blue) of the symmetric mode in a pair of waveg- 
uides in the backwards wave regime, with a negative overlap. 



erally constitutes the largest contribution to the overlap 
integral, especially close to the cutoff frequency. Its pro- 
file E x (z) is sketched in Fig. p] (Supplemental) (a). The 
E x field of an antisymmetric mode is itself symmetric, 
and it is continuous at the metal-dielectric interfaces. As 
a result, j \E X (z) E x (z — d) dz > 0, just like the overlap 
integral of E Z1 and therefore we know that n must be 
positive. 

Finally, the E x field profile of the symmetric mode in a 
pair of waveguides (i.e., the odd superposition of the an- 
tisymmetric modes of the individual MDM waveguides) 
is depicted in Fig. [6] (Supplemental) (b). Since the sym- 
metry is opposite to that of the backwards wave in a 
single MDM waveguide, in this case k < 0. This now 
fully explains why refraction of both phase and energy 
are negative in a metal-dielectric multilayer stack with 
alternating thick and thin metal films. 

The above analysis has again illustrated how mode 
symmetry governs the sign of the coupling constant. 
In particular, in metal-dielectric waveguides (where the 
transverse component E z switches sign at the core- 
cladding interface), symmetric modes yield n < and 
antisymmetric modes yield n > 0. We note that in all- 
dielectric waveguides — in which E z does not switch sign 
across an interface — the relationship between symmetry 
and the sign of the coupling constant will in most cases be 
reversed. The connection is however less general in that 
case, since the transverse and longitudinal components 
have opposing contributions to the overlap integral. 



Isofrequency contours in the presence of absorption 

It is important to evaluate the effect of realistic absorp- 
tion losses and compare the results to the lossless case. 
We use e m = —4.43 + 0.2JJ in the example described in 
the main text at Ao = 400 nm [5 , to plot the wavevec- 
tor diagrams in Fig. [7Fa) for two cases; taking k z real 
and k x complex (blue curve, relevant in the context of 
propagation and damping along the layers, and refrac- 
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FIG. 7: (Supplemental). Calculated wavevector diagrams for 
the example structure described in the main text, including 
absorption in the Ag, taking either k x or k z to be complex. 
(a) Re (fee) as a function of real k z (blue) and Re(k z ) as a 
function of real k x (red), (b) Im (k x ) as a function of real k z 
(blue) and Im(k z ) as a function of real k x (red). The grey 
dotted curves represent k x + k 2 z — k$ . 



tion at an interface of the metamaterial parallel to the 
7/z-plane) and taking k x real and k z complex (red curve, 
useful in the context of propagation normal to the layers 
and refraction at an interface parallel to the x?/-plane). 
Because the red curve still follows the n — — 1 contour 
(the grey dotted curve) at small \k x \ and the blue curve 
does the same at small \k z \, it can be concluded that 
losses do not significantly affect the value of the effective 
index, nor its isotropic character. In Fig. [7] (Supplemen- 
tal) (b), the imaginary part of the wavevectors is plotted 
for both cases. It shows that Im(/c 2 ) is relatively small 
for — ko < k x < &o, where the mode has a propagat- 
ing character. Outside this regime, the mode is largely 
evanescent. Likewise, Im (k x ) is small for —ko<k z <ko, 
and has a maximum value at k z = ir/d. The figure of 
merit FOM = Re (|k|) /Im (|k|) is only unambiguously 
defined along the x and z axes, since in those directions 
the symmetry of the geometry dictates that k and S are 
entirely parallel or antiparallel. The FOM varies from 
14.7 along the x direction (obtained by taking k z = 
and k x complex) to 8.7 along z (with k x = and k z 
complex) . 



Calculation of coupling efficiency 

The coupling efficiency of a wave incident from air 
on an interface of the metamaterial parallel to the yz 
plane can be evaluated in a finite- difference time-domain 
(FDTD) simulation by integrating the power flow inside 
the metamaterial at a depth of 1 /im below the interface. 
At this depth, only the propagating backwards mode ex- 
ists in the waveguide. Taking into account the known 
propagation losses over the first 1 /im, the fraction of 
energy coupled to the mode of interest is calculated. 

By varying the layer thicknesses in the array of coupled 
waveguide pairs, in particular the relative thicknesses of 



the two metal films, the coupling efficiency can be con- 
trolled to some extent. For example, at a free- space 
wavelength of 474 nm, a coupling efficiency of 19.6% was 
found for a stack consisting of Ag (using permittivity 
from ref. [5J and a dielectric with a refractive index of 
3.12. The dielectric thickness was 44 nm, and the thick- 
nesses of the alternating metal films were 24 and 90 nm, 
respectively, resulting in an effective index of —2.1 at 
normal incidence. We note that we limited the investi- 
gated structures to those consisting of only one type of 
metal and one type of dielectric, and we have chosen the 
thickness of the two dielectric films in the four-layer unit 
cell to be equal. Including different materials and allow- 
ing the thicknesses of the two dielectric layers to be dif- 
ferent would introduce additional degrees of freedom to 
further optimize impedance matching. Moreover, it may 
be possible to improve impedance matching by properly 



structuring the interface of the metamaterial [6 . 
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